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Abstract 


In  the  thermal  generatim  of  nuclear  magnetic  resonance  phenomena,  such 
as  is  found  in  flames  or  hot  eidiausts.  energy  is  lost  during  the  transition  from 
a  hi^er-  to  a  lower-energy  nuclear  mapietic  quantum  state.  The  interaction 
that  takes  place  is  explained  in  simple  theoretical  terms,  and  the  effects  of  re¬ 
sonance  broadening  by  spin-spin  interaction  and  an  inhomogeneous  magnetic 
field  are  cmisidered.  The  amount  of  energy  lost  during  the  cooling  of  a  hot  jet 
exhaust  was  calculated  using  the  environmental  ma^ietic  field  of  the  earth. 

On  the  assumption  that  the  energy  was  radiated  by  means  of  electromagnetic 
waves,  attempts  were  made  to  detect  this  energy  loss. 
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NUCLEAR  A^GNETIC  RESONANCE  IN  FLAMES 


Introduction 


The  phenomenon  of  nuclear  magnetic  resonance  is  probably  best  recognized 
in  the  fields  of  organic  chemistry  and  nuclear  physics.  In  the  years  following 
its  discovery  by  Purcell  and  Bloch  in  1946.  the  phenomenon  has  been  applied 
in  the  chemical  analysis  of  solids,  liquids,  and  gases.  Comparatively  little 
consideration,  however,  has  been  given  to  the  thermal  generation  of  the  reson¬ 
ance  which  should  take  place  as  a  gas  undergoes  a  rapid  change  in  temperature. 
Any  rapid  change  in  temperature,  whether  it  is  an  increase  or  a  decrease, 
should  induce  resonance  just  as  immersion  of  the  medium  in  an  oscillatory  mag¬ 
netic  field  induces  resonance. 

A  simple  explanation  of  nuclear  magnetic  resonance  will  provide  a  necess¬ 
ary  background,  and  the  theory  will  be  expanded  to  apply  to  some  medium  under¬ 
going  a  rapid  decrease  in  temperature.  With  this  background,  the  amount  of 
energy  transfer  possible  for  several  commercial  jet  engines  using  JP-4  pro¬ 
pellant  can  be  determined  quantitatively. 

Attempts  were  made  to  detect  any  electromagnetic  radiation  produced 
during  the  cooling  of  the  hot  jet  exhaust  gas.  and  the  method  and  results  are 
described  in  detail.  A  simple  approach  to  the  thermal  generation  of  nuclear 
magnetic  resonance  is  presented,  and  the  problems  encountered  in  practical 
application  are  discussed. 


1.  Th«or«tical  Aspects 

The  property  of  importance  in  nuclear  magnetic  resonance  is  the  nuclear 
spin  I,  or  the  circulation  of  the  nuclear  charge.  Such  charge  circulation  gen¬ 
erates  a  magnetic  field;  hence,  the  nuclear  magnetic  moment  vector  p  arises. 

The  nuclear  magnetic  moment  vectors  are  quantized,  or  can  take  on  only 
certain  specified  average  values.  Suppose,  for  example,  that  the  nuclear  spin 
I  equals  1/2.  (I  ■  1/2  will  be  used  throughout. )  The  two  possible  orientations 
of  the  magnetic  moment  n  when  it  is  immersed  in  a  magnetic  field  are 
shown  in  Fig.  1.  The  two  possible  quantum  states  are  m  ■  t-l/2  and  m  *  -1/2. 
where  -I<  m  <  I. 


Figure  1. 
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Two  Possible  Orlentstioas  of  The  Magnetic  Monent  m 
When  it  is  Immersed  In  a  Magnetic  Field 
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In  the  absence  of  a  magnetic  field,  there  is  no  particular  preference  for  one 
state  over  the  other,  but  in  the  case  of  an  impressed  magnetic  field,  the  ■•-1/2 
state  becomes  the  more  favorable.  This  is  simply  shown  by  likening  the  mag¬ 
netic  moment  vector  to  a  compass  needle.  Orientation  in  the  same  direction  as 
the  impressed  field,  m  =  -*-1/2,  corresponds  to  a  normal  compass  needle  orien¬ 
tation.  Orientation  in  the  opposite  direction,  m  >  -1/2,  corresponds  to  the 
compass  needle  rotated  180*.  The  latter  state  requires  work  to  be  expended; 
thus  the  -1/2  state  is  the  more  energetic,  hence  the  less  favorable,  state  for 
occupation.  If  it  were  not  for  thermal  agitation,  the  tendency  to  assume  the 
more  favorable  +1/2  state  would  be  unopposed. 


If  the  population  of  the  lower  energy  state  is  to  be  opposed  by  thermal  agi¬ 
tation,  there  must  be  some  mechaiiism  for  the  transfer  or  coupling  of  energy 
from  the  surrounding  medium  (the  lattice)  into  the  lower  energy  spin  quantum 
state.  This  mechanism  is  provided  through  the  precession  of  the  magnetic 
moment  about  the  steady  impressed  magnetic  field.  If  there  is  to  be  a  transfer 
of  energy,  an  oscillatory  magnetic  field  of  the  same  phase  and  frequency  as 
that  of  precession  must  be  impressed  upon  the  medium.  This  condition  produces 
a  magnetic  torque  which  tends  to  ' flip'  or  rotate  the  magnetic  moment  vector 
into  the  higher  energy  orientation,  as  in  the  case  for  a  simple  compass.  Thus, 
energy  is  absorbed  from  the  source  that  produces  the  oscillatory  magnetic  field. 


To  find  the  population  distribution  of  energy  states  under  the  condition  of 
thermal  equilibrium  for  the  spin  temperature  and  the  temperature  of  the  lattice, 
the  Boltzmann  distribution  law  is  applied  with  the  nuclear  magnetic  moment  n, 
the  steady  magnetic  field  strength  H^,  and  the  absolute  temperature  T; 
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where  is  the  number  of  nuclei  in  the  given  m  state,  is  the  total  number 
of  nuclei,  ia  the  energy  associated  with  the  given  m  state,  and  k  is  Boltz¬ 
mann's  constant. 
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where  is  now  written 

Mh  '  >1  ™ 

The  nuclear  gyromagnetic  ratio  is  and  is  Planck's  constant . 
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In  the  case  of  different  temperaturea,  aay,  along  a  flame 

from  a  higher  to  a  lower  temperature  (T^  >  Tj ) : 

For  T,  Vl  ^  **o 

N  ah*. 

T,  '  ^  • 

VlMp 

N  2kT, 

For  Tj  T  -  -ar  • 


(N^l/a  »  Tj  ^  ^**+1/2*  Tj  • 

Thua  for  a  change  of  temperature  as  described  here,  the  number  of  nuclei  in 
the  *1/2  state  is  increased.  Correspondingly,  energy  must  be  lost  to  the  sur¬ 
roundings  because  there  has  been  a  decay,  as  it  were,  from  the  -1/2  to  the 
-*■1/2  state.  The  number  of  nuclei  that  have  experienced  this  decay  is  indicated 

by  the  increased  p<^lation  of  the  lower  energy  state.  .  . 

Yl  Ji  Yii\  ^ 

/  TVT'  TLT .  1 

‘"■'".w'tj-'Wt,-  ^ 

Expanding  exponentially. 


The  rate  of  energy  loss  is  proportional  to 
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In  the  steady  state  g  ^  =  0  .  Then 

=  i_I_y  +_LTv  +  -LI  V 

dt  8x*  By  y  8z* 
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In  the  case  of  a  flame,  it  is  necessary  to  know  the  velocity  components  and 
temperature  gradients,  along  with  and  N^,  to  determine  the  amount  of  energy 
loss. 

In  the  case  of  an  exhaust,  such  as  a  Jet  engine  exhaust,  and  considering 
only  one  -dimensional  (axial)  flow: 


where  and  are  the  exit  and  ambient  temperatures  respectively,  and  v^ 
is  the  average  velocity.  Three  specific  cases  of  interest  are  enumerated  in 
Appendix  A.  Figure  2  shows  the  relationship  between  power  loss  in  watts  to  the 
temperature  gradient  for  these  three  cases. 

Various  authors  (including  this  one)  have  assumed  that  the  power  lost  is 
immediately  radiated  to  the  lattice  by  means  of  electromagnetic  radiation  of 
the  proper  frequency.  This  energy  may  simply  be  returned  to  the  lattice  by  the 
reverse  of  the  coupling  mechanism  discussed  earlier.  The  energy  lost  would 
then  be  realized  as  a  heating  of  the  lattice. 

The  preceding  discussion  was  concerned  with  a  simple  approach  to  a  very 
complex  interaction.  There  are  many  other  factors  that  could  have  grave  con¬ 
sequences  so  far  as  the  production,  loss,  and  detection  of  magnetic  resonance 
energy  in  flames  are  concerned. 

The  foremost  factor  is.  of  course,  that  the  magnetic  spin  energy  loss  may 
be  simply  realized  in  the  form  of  thermal  energy  without  any  subsequent  electro¬ 
magnetic  radiation  at  all.  The  extent  of  the  spin-lattice  interaction  must 
therefore  be  considered  more  in  detail.  The  spin  is  the  spin  of  the  nuclei  and 
the  lattice  is  the  framework  or  material,  whether  gas,  liquid,  or  solid,  within 
which  the  nuclei  are  found.  Present  opinion  is  that  while  this  interaction  is  so 


Figure  2.  Power  Lose  vs  Tempersture  Orsdient  for  Selected  Jet  Engines 
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small  that  spin  temperature  and  lattice  temperature  may  be  considered  inde¬ 
pendently,  what  interaction  there  is  between  the  two  systems  seems  to  bring 
them  into  equilibrium  at  the  same  temperature. 

Second,  equilibrium  statistics  are  being  applied  to  a  system  that  most 
assuredly  is  not  in  equilibrium.  Since  the  proper  statistical  representation  is 
not  known,  the  results  are  q;>en  to  question. 

Third,  the  resonance  broadening  of  spin-spin  interaction  comes  about  be¬ 
cause  many  of  the  nuclei  have  magnetic  dipole  moments,  hence  magnetic  fields, 
associated  with  them.  The  effect  of  these  magnetic  fields  is  to  place  the  nuclei, 
upon  which  the  resonance  phenomenon  is  dependent,  in  a  very  heterogeneous 
state:  one  nucleus  'feels'  a  magnetic  field  of  a  certain  strength,  while  a  neighbor 
'feels'  a  different  strength,  ana  so  on.  Thus,  a  band  of  frequencies  is  produced 
instead  of  a  single  frequency.  Normal  resonance  techniques  could  be  applied 
to  a  sample  in  equilibrium  and  the  extent  of  this  effect  determined  to  an  order 
of  magnitude,  provided  it  is  within  the  resolving  power  ot  the  spectrometer. 

A  theoretical  approximation  shows  this  effect  to  be  negligible  in  the  gas  under 
consideration.  (See  Appendix  B. ) 

A  fourth  factor  is  the  homogeneity  of  the  magnetic  field  in  which  the  entire 
system  is  immersed.  This  need  not  be  discussed  further  than  to  indicate  that 
a  change  ot  1  cps  in  the  resonant  frequency  for  protons  is  produced  by  a  change 
in  magnetic  field  of  approximately  2&  gamma  (1  gamma  *  10*^  gauss>  .  The 
same  broadening  effect  as  for  spin-spin  interaction  would  be  noted. 

Fifth,  the  radiation,  if  any,  would  be  incoherent. 

Sixth,  no  consideration  is  given  to  tne  time  of  relaxation  from  the  higher  to 
the  lower  energy  state.  It  is  assumed  that  the  relaxation  takes  place  at  the 
same  rate  as  the  thermal  relaxation  dT/dt  =  v^  (dT/dx)  .  This  may  not  be  the 

case. 

2.  ExperiiTMntal  Report  and  RmuRs 

Based  on  the  precedmg  calculations,  several  attempts  were  made  to  detect 
nuclear  magnetic  resonance  as  supposedly  produced  in  Jet  exhausts.  The  cri¬ 
tical  assumption  was  that  the  energy  associated  with  the  resonance  phenomenon 
would  be  radiated. 

Measurements  were  made  under  line-of-sight  conditions  from  the  Katahdin 
Hill  field  site  of  the  Astrosurveillance  Sciences  Laboratory  to  the  jet  engine 
static  test  stand  (approximately  one-half  mile  away)  .  The  frequency  of  reson¬ 
ance  could  be  only  approximated  and  the  bandwidth  was  unknown.  For  these 
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reasons,  a  sweep  frequency  technique  was  used  in  order  not  to  lose  any  infor¬ 
mation  nor  to  have  long,  noisy  engine  runs.  Furthermore,  the  expected  signal 
had  a  low  intensity  level,  but  the  new  electronic  analog  multiplier  developed  by 
the  Electronic  Development  Section  of  the  Astrosurveillance  Sciences  Labor¬ 
atory  afforded  an  extra  gain  characteristic  in  signal  detection.  A  block  diagram 
of  the  entire  syetem  is  shown  in  Fig.  3  and  an  equipment  list  and  specifications 
are  found  in  Appendix  C. 


The  process  involved  in  the  system  was  pseudo  cross  correlation: 


Let  Acoswt  be  the  signal  input  from  the  sweep  oscillator  where  for  the 
moment  uf  f(t)  and  A  is  never  a  function  of  time.  Let  A^cosw^t  be  the  signal 
to  be  detected  where,  again  for  the  moment.  A^  and  j  f(t)  .  The  multiplier 
multiplies  and  integrates  the  two  signals  over  a  given  time.  The  mathematics 
are: 


( 1 )  Multiplication 

A  cos  (lit  A^  cos  oi^t  >  AA^  cos  ait  cos  b).  t 

(2)  Integration 

/T  [sin  )t 

cos  «t  cos  dt  •  AAJ-j^— I  ■ 

i 

AAj  r  sin  <  w  -  u)^  jr  sin  (w  ^  )  T  1 


sin  (w+Wj  It]"^ 

*  Jo 


This  result  may  be  analyzed  by  saying  that  as  w  the  first  term  goes  to 
infinity.  The  second  term  is  always  very  small  in  comparison  and  may  be  neg¬ 
lected.  Therefore, 

AA.  sin  2x(f-f,)  T 
A(T)  =  - - - 1 - 

2x  <f-fj ) 


In  the  actual  experiment  the  following  parameters  were  used  for  a  magnetic 
field  strength  of  0.  56  gauss: 


f  «  r  +  k  T 
o 

{  ranges  from  2290  to  2560  cps  as  r  ranges  from  0  to  30  sec 
f^  >  2290  cps 
2560  -  2290  30k 

f  •  2290  9  T 

/AA,\  sin  2x  (2290 9t  -  f.)  T 

A(T)  -( - i) - i - 

V  2  /  2w  (2290  +  9t  -fj) 
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As  an  example  of  the  theoretical  response,  *  2400  cps. 


A(t) 


.  sin  2 ir(  9  r  -  110) 

Sin  2w  (9t  -110)  T 

'  2ir  (9t  -110) 

^4  ) 

9  T-  no 

This  response  was  filtered  to  pass  only  those  frequencies  balow  10  cps.  At 
that  portion  of  the  equation  near  f^  ■  2400  cps  covered  by  T  ranging  from  18  to 
22  sec,  the  response  would  be  that  shown  in  Fig.  4.  Figure  5  is  a  trace  of  an 
actual  pickup  by  the  system. 

The  disadvantage  of  this  system  was  that  the  frequencies  had  to  be  read 
from  a  calibrated  scale.  This  was  tedious  since  the  speed  of  the  sweep- 
oscillator  varied  from  one  sweep  to  the  next,  requiring  the  construction  of 
numerous  scales.  Furthermore,  human  Judgment  was  needed  in  the  case  of 
small  signal  possibUities. 

Data  read  off  was  recorded  as  a  panoramic  frequency  plot  showing  number 
of  occurrences  vs  frequency  for  times  when  the  engine  was  in  operation  and 
times  when  it  was  not.  Figures  6,  7.  8,  and  9  show  the  results. 

The  only  indications  of  possibly  abnormal  activity  were  on  the  KC-135  test. 
Two  bands  of  activity  were  observed,  one  at  about  2360  to  2365  cps  and  the 
other  at  about  2475  to  2485  cps.  Previous  noise  measurements  showed  expected 
peaks  at  2320,  2385,  2440,  2485  (weak)  ,  and  2560  cps.  These  noise  measure¬ 
ments  were  verified  because  each  of  these  frequencies,  with  the  exception  of 
2495  cps.  was  usually  detected.  A  further  problem  was  that  the  nuclei  could 
not  be  controlled  so  that  radiation  would  occur  in  one  of  the  observed  quiet  bands. 


Attempts  were  made  to  study  the  amplitude  of  received  signals  from  the 
characteristic  of  the  system  output.  Trends  were  found  that  indicated  the  pos- 
sibilityof  calibration  by  means  of  the  central  maximum  deflection  and  to  a  lesser 
extent  the  Vidth  of  the  observed  characteristic.  The  equipment  would  have 
needed  greater  refinement,  however  for  more  precise  calibration. 


The  possibility  of  a  high  incidence  of  error  in  Judgment  made  mandatory 
machine  analysis  of  the  data.  The  final  tests  were  made  with  the  system  shown 
in  Fig.  10.  The  operation  of  this  system  is  simple.  The  received  signal  is 
amplified  and  fed  into  a  mixer  where  it  is  beat  against  a  local  oscillator  set  at 
ai^roximately  the  center  frequency  of  the  band  of  interest  (2300  to  2550  cps)  . 
The  difference  is  then  passed  by  the  low  pass  filter,  amplified,  and  recorded 
at  slow  speed  on  magnetic  tape.  The  tape  playback  is  speeded  up  (in  this  case 
four  times  to  fit  the  Vibralyzer  range  of  5  to  500  cps)  .  fed  into  the  Vibralyzer 
and  subsequently  analyzed.  No  peculiarities  were  noted,  and  the  results  were 
negative. 
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3.  Conclusions 

Although  there  is  probably  a  reaonance  energy  loss  due  to  a  change  of 
nuclear  magnetic  energy  state  from  the  higher  level  to  the  lower  level,  it  was 
not  possible  to  detect  the  loss  by  means  of  radiation  techniques.  This  may  be 
attributable  to  the  inhomogeneous  conditions  (magnetic  field)  under  which  the 
measurements  were  made.  In  only  one  case  out  of  seven  test  runs  was  there 
the  slightest  possibility  of  detection  of  any  sort  of  Irregularity,  and  that  occurred 
in  two  bands. 

Under  conditions  of  strictly  controlled  magnetic  field  and  thermal  gradient 
with  a  low  noise  background,  the  energy  loss,  if  in  the  form  of  electromagnetic 
radiation,  might  be  detected.  In  fact,  experiments  under  highly  controlled  con¬ 
ditions  should  answer  the  several  problems  posed  earlier. 
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Fifur*  4.  TiMorctlcal  System  Output 


MMCM 


MWCW 


I 


13 


ENGINE  TYKt  J-33 


TEST  NO.  1  TEST  NO.  2 


Figure  6,  Number  of  Occurrences  vs  Frequency,  21  November  1960 
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APPENDIX  A 

Engine  Characteristics 


Engine  type:  J-33 

Thrust:  5000  lb 
Fuel  type:  JP-4 
Hydrogen  content  of  fUel:  13.  6% 
Fuel  rate:  1 . 44  lb/  sec 


Proton  rate:  Fuel  rate  X  Hydrogen  content  of  fuel  „ 
“WelgBr  of  a  single  hydrogen  nucleus  *  ^o 


N 


1.44  Ib/sec  X  453.8  gm/lbX  0.136  _ _ 25 

=  — — ' - .24  ' - *  5.  32X  lO*"*  proton/sec 

1.  67  X  10  gm/nucleus 


Exhaust  temp:  -  1240*F  ■  940‘K 

Ambient  temp:  »  300*K 


Exhaust  velocity  (axial)  : 


y  ,  Thrust  X  yceleratlon  of  gravity 
E  Fuel  rate 


,  5000  lb  X  32. 2  ft/sec* 

- 1.44  1b/sec  ' - 

Ambient  velocity:  v^  »  4. 58  X  10 

Proton  gyromagnetic  ratio:  Yi  * 
Population  shift: 


■=  11. 1  X  10*  ft/sec  * 
*  cm/sec 

2.  675  X  10*/sec/gauss 


338  X  10* 


cm/sec 


dN^/(5.32X10^^\ 
dt  *  \  4  / 

(2.  675X  10*/sec/gauss)(l.  05X  10~^%rg/sec)(0.  56gatmsXl67X10*cm/sec)  /dT\ 
(1.  38  X  10*l*erg/deg>  (940  X  300  deg*)  'dx  ' 


dN 

W 


8.97  X  10^®(-I^)/  sec 
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Energ>/ shift:  hi/  ■  1.  57  X  10~^^  ergs 
Power  loss:  P  •  1. 42  X  10'^^  ^^^wstts 


Engine  type:  J-57 


Thrust:  10.000  1b 
Fuel  type:  JP-4 

Hydrogen  content  of  Aiel:  13.0% 

Fuel  rate;  3.  6  Ib/sec 
Proton  rate: 

N  .  3.6  1b/secX  453gm/lbX  0.136  .  ^  33  ^  proton/sec 

1.  67  X  10'*  gm/nucleus 

Exhaust  temp:  -  1000*F  ■  813*K 

Ambient  temp:  ■  300*K 


Exhaust  velocity  (axial)  : 

„  .  10,  000  lb  X  32. 2  ft/sec* 
'^E  - 3.  B1b7sic'  .  - 


8.6  X  lO^ft/sec 


274  X  lO^cm/sec 


Ambient  velocity:  v^  >  4. 58  X  10^  cm/sec 

Proton  gyromagnetic  ratio:  Vj  ■  2.675  X  10^/sec/gauss 

Population  shift: 

dN  /(I.  33X10^*)  y2.  675Xi0*)(l.  05X10~^'^)(0.  56)(135XlO*)  /^\ 
dt  \  4  /  (1  38X10"^*)(813XS00)  Vdx  / 


gi.2.10X10'»(^)/sec 

Energy/shift:  hi/«  1.57  X  10*^^  er|^ 

Power  lose:  P-  3. 30  X  watts 

Engine  type:  Four  J-S7's  (KC-135) 

Characteristics  all  4X  J-57's 

Power  loss  KC-135:  P-  1. 32  X  lO"'®  watts 
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APPENDIX  B 


A  theoretical  approximation  of  the  spin- spin  interaction  is  as  follows: 

Assume  the  ideal  engine  in  which  the  exhaust  pressure  equals  the  ambient 
pressure.  The  mean  free  path  in  the  exhaust,  based  on  hard  sphere  considera¬ 
tions  and  a  Maxwellian  velocity  distribution,  is 

KT 

V2  pxd^ 

where  p  is  the  pressure  and  d  is  the  diameter  of  the  molecules.  Combustion 
products  are  principally  COj  and  H^O.  Taking,  then,  the  length  of  the  COj 
bond  as  d, 

d  *  1.42  X  10  '^cm 

d^  *  2.02  X  10'^®cm^ 

p  *  1.  01  X  10*dynes/cm^ 

T  •  lOOOT  «  813*K 

,  _ (1. 38X10' ^^dyne  cm/deg) (813  deg) 

(1 . 414)(  1. 01 X  I0^dyne/cm^> (3. 14) (2.  02 X 10' ^®cm* 

X  »  1.24  X  10"®cm 


The  local  field  produced  by  a  nucleus  having  the  magnetic  moment  of  a  proton 
/ip  *  1.41  X  1  o' ^^erg/ gauss 


is  given  as 
“local 


1.41  X 

(1.24X10'^  ^ 


7.34  X  10*'^  gauss 


The  local  field  is  negligible. 
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APPENDIX  C 

System  Characteristics 


SyBtem  No.  1 
Antenna 

3  ft  by  4  ft.  274  turns.  No.  18  wire 
Bandwidth  at  2400  epa  «  200  cps 

Transistorized  amplifier 
Gain:  SO,  000 

Bandpass:  SOO  cps  to  6  keps 
Flat  from  1 . 5  keps  to  3  keps 
Equivalent  rms  input  noise  «  0.  S 
Saturation  «  30  fiv 

Multiplier 

See  Frederick  Slack.  Electronic  Analog  Multiplier.  AFCRL  541,  July  1961 
Recorder 

Sanborn  Twin-Viso  Model  60-1300 

Sweep  frequency  oscillator 

Modified  General  Radio  Type  1304B 

Receivers 

Stoddardt  NM-40A  Radio  Interference  Field  Intensity  Meter 

System  No.  2 

FM  Recorder 

Ampex  FR-IOOA 

Vibralyzer 

Kay  Lab  "Vibralyzer" 
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